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ABSTRACT 

Excited-state OH maser emission has previously been reported in the circumstellar envelopes of only two 
evolved stars: the Mira star AU Geminorum and the hypergiant NML Cygni. We present Very Large Array 
(VLA) observations of the 1665, 1667, and excited-state 4750 MHz mainline OH transitions in AU Gem 
and Expanded Very Large Array (EVLA) observations of the excited-state 6030 and 6035 MHz OH mainline 
transitions in NML Cyg. We detect masers in both mainline transitions in AU Gem but no excited-state emission 
in either star. We conclude that the excited-state OH emission in AU Gem is either a transient phenomenon 
(such as for NML Cyg outlined below), or possibly an artifact in the data, and that the excited state OH emission 
in NML Cyg was generated by an episode of enhanced shock between the stellar mass-loss and an outflow of 
the Cyg OB2 association. With these single exceptions, it therefore appears that excited-state OH emission 
indeed should not be predicted nor observable in evolved stars as part of their normal structure or evolution. 
Subject headings: masers — circumstellar matter — stars: individual (AU Gem, NML Cyg) — stars: late type 
— stars: mass loss 
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1. EVITRODUCTION 

The topic of hydroxyl (OH) maser pumping in the circum- 
stellar envelopes of evolved stars has received increased fo- 
cus over recent years due to observational and theoretical ad- 
vances. On the observational side, detections of far infrared 
OH lines toward evolved stars have confirmed the role of ra- 
diative pumping routes and allowed for estimation of pump 
efficiencies (Sylvester et al. 1997; Neufeld et al. 1999; He & 
Chen 2004; He et al. 2005). On the theoretical side, the re- 
cent Gray et al. (2005) model demonstrates the complexity of 
OH pumping routes and is quite successful at explaining ob- 
served OH maser properties from OH/IR stars. Pump models 
generally do not predict detectable excited-state emission in 
circumstellar envelopes of evolved stars, either due to a lack 
of inversion or insufficient optical depth (e.g., Elitzur et al. 
1976; Bujarrabal et al. 1980; M. D. Gray 2006, private com- 
munication). For a comprehensive overview on circumstellar 
shells around evolved stars we refer to Habing (1996). 

Excited-state OH masers have been sought in the 4.8 and 
6.0 GHz transitions toward the circumstellar environments of 
a variety of evolved high mass losing stellar sources by many 
authors (Thacker et al. 1970; Zuckerman et al. 1972; Baudry 
1974; Rickard et al. 1975; Claussen & Fix 1981; Jewell et al. 
1985; Desmurs et al. 2002; Fish et al. 2006). These searches 
universally failed to detect excited-state emission except in 
two circumstellar environments. Claussen & Fix (1981) re- 
port on a 5(7 detection of a 4750 MHz maser in AU Gem, 
a Mira variable. Subsequent observations by Jewell et al. 
(1985) failed to confirm this emission, but the RMS noise 
level of their observations does not rule out the 100 mJy de- 
tection of Claussen & Fix (1981). 

The second evolved star with an excited-state maser detec- 
tion is NML Cyg, a supergiant and recently more often re- 
ferred to as a hypergiant (van Genderen et al. 1982; Schus- 
ter et al. 2006). Zuckerman et al. (1972) report a clear (^ 

' National Radio Astronomy Observatory, 1003 Lopezville Rd., Socorro, 
NM 87801, lsjouwer@nrao.edu, vfish@nrao.edu, mclausse@nrao.edu. 
^ Jansky Fellow. 

Department of Physics and Astronomy, University of New Mexico, 800 
Yale Boulevard NE, Albuquerque, NM 87131, ylva@unm.edu. 



10 cr) detection at 6035 MHz and possible emission ('^ 3 a) 
at 6030 MHz as well. Their spectra indicate strong (2.2 Jy) 
6035 MHz maser emission at Vlsr « 4 kms"'. This is close 
to the stellar velocity of about (± 4) km s"' (Zhou 1981; 
Diamond et al. 1984), usually found as the average velocity 
of the outer boundaries or emission peaks of the characteris- 
tic double -peaked 1612 MHz OH emission or as the center of 
the SiO emission in OH/IR stars (Bowers et al. 1983; Habing 
1996). With ample sensitivity to follow up on this initial de- 
tection, Jewell et al. (1985) did not detect any emission from 
NML Cyg a decade later In 1999, Desmurs et al. (2002) also 
did not detect the emission at 4 km s~' but note possible weak 
(^ 20 mJy) emission at -17 kms~\ corresponding in veloc- 
ity to an inner 1612 MHz peak near -18 kms"' (e.g., Herbig 
1974; Engels 1979). The original spectra of Zuckerman et 
al. (1972) are suggestive of weak emission at both 6030 and 
6035 MHz near this velocity, but not at significant levels com- 
pared to their RMS noise. 

The confirmed detection of excited-state OH masers in the 
circumstellar environments or ejecta of evolved stars would 
provide a challenge to modern pumping models, or at least 
highlight novel portions of pumping phase space heretofore 
not considered. To date, neither report of a detection of an 
excited-state OH maser in an evolved star has been confirmed, 
which in view of the pumping models therefore needs to be 
followed up. In this Letter, we present new observations of 
the OH masers in AU Gem and NML Cyg in the excited-state 
OH lines originally reported as detections, as well as ground- 
state 1665 and 1667 MHz mainline OH emission in AU Gem. 

2. OBSERVATIONS 

AU Gem and NML Cyg were observed with the Very Large 
Array (VLA) using the settings as outlined in Table 1 . The AU 
Gem observations occurred in DnC configuration. Due to the 
Expanded VLA (EVLA) upgrade (McKinnon & Perley 2001; 
Ulvestad et al. 2007), only 22 antennas were available, includ- 
ing three EVLA antennas. All these antennas were used. 

For NML Cyg we were able to profit from the special call 
for proposals in 2007 April for using the new 5 GHz (C-band) 
receivers on the EVLA. These receivers are part of the EVLA 
upgrade and allow observing at a much wider frequency 
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TABLE 1 
Observational summary table 

Observing Date Line Au AV Beam RMS 

(MHz) (kHz) (kms-') ("x",°) t 



AU Gem, centered at Vlsr ~ 10 km s ' 



2007 Jan 25 


1665 


195 


0.28 


40.3x21.5,73 


2007 Jan 25 


1667 


195 


0.28 


40.4x20.4, 72 


2007 Jan 25 


4750 


781 


0.19 


13.0x12.4,94 


NML Cyg, centered at Vlsr 


« 0km s 


-1 




2007 May 26/30 


6030 


1563 


0.61 


0.56x0.23, 118 


2007 May 26/30 


6035 


1563 


0.61 


0.58x0.23, 118 




10 20 10 20 

LSR Velocity (km s"') 



Fig. 1. — Left: Spectra of the 1665 MHz emission in AU Gem. RCP 
emission is shown in bold and LCP in normal weight. The Stokes I spectrum 
has been shifted vertically for clarity. Grey dots indicate channels in which 
emission is detected at the phase center at more than 5 times the RMS noise. 
Right: Spectra of the 1667 MHz emission in AU Gem. 

emission is stronger in each transition than left circular polar- 
ized (LCP) emission. The second peak near 15 kms"' corre- 
sponds to the broader redshifted shoulder of the emission in 
the lower-resolution Nguyen-Q-Rieu et al. (1979) spectrum. 
We detect weaker features at higher LSR velocity as well, in- 
cluding emission in the edge spectral channels at ~ 28 km s~' 
at 1667 MHz, but we do not detect the weaker (^ 70 mJy) 
emission at km s"' claimed by Nguyen-Q-Rieu et al. (1979). 
If we assume that the 1665 and 1667 MHz mainline OH-peaks 
at the systemic velocity of 13 kms~' define the stellar veloc- 
ity, and if we conclude that the Nguyen-Q-Rieu et al. (1979) 
~ km s"' 1667 MHz weak feature is convincing, then these 
features could indicate an expanding OH shell with an expan- 
sion velocity of about 14-15 kms"'. If we were to dismiss 
the Nguyen-Q-Rieu et al. (1979) ~ kms"' 1667 MHz fea- 
ture, the systemic and expansion velocities would be 20 and 
8 kms"', respectively. Both expansion velocities are reason- 
able for Mira-type stars with the latter for a less optically thick 
(lower metallicity) shell (Habing 1996), but the uncertainties 
in the spectra do not allow a firm conclusion on this. In this re- 
spect it is unfortunate that AU Gem has never been detected in 
the 1612 MHz OH transition (Fix & Weisberg 1978; Nguyen- 
Q-Rieu et al. 1979; Olnon et al. 1980) nor shows any H2O nor 
SiO emission (Nyman et al. 1986). All emission is consistent 
with being point- like at the 40" x 20" resolution of the VLA in 
DnC configuration, which is as expected since a typical Mira 
OH shell size (10'^ cm, from Herman & Habing 1985) would 
subtend an angle of less than 0''3 at the distance of AU Gem 
(2.4 kpc, from Nguyen-Q-Rieu et al. 1979). 

As to the 4750 MHz transition, we do not detect any emis- 
sion at a 3 (7 noise level of 10 mJy beam"' in AU Gem. 

3.2. NML Cyg 

No emission was found over 4.2tT (34.2 mJy beam"') in 
the 6030 MHz data, nor was any emission found over A.9a 
(35.3 mJy beam"') in the 6035 MHz data within the beam 
of the Jewell et al. (1985) observations. We cannot rule out 
tentative (>4.8(t single channel) 6035 MHz line emission at 
-21.2, +16.4 and +20.6 km s"' (Fig. 2). This would support 
the notion of tentative 6035 MHz OH features at » -17 km 
s"' (Desmurs et al. 2002; Zuckerman et al. 1972) and could 
be related to a new episode of high mass loss or compression 
of circumstellar material close to the star as traced by the col- 
lisionally excited SiO maser (Boboltz & Marvel 2000). How- 
ever, we would not claim such emission without extensive 
high-sensitivity follow-up observations, e.g., when all VLA 
antennas have new C-band receivers installed (in 2010). 



t Typical RMS noise in a spectral channel in mJy beam 
"Assuming a flux density of 2.0 Jy at 1.7 GHz for J074 1+3 12 
Assuming a flux density of 1.6 Jy at 4.8 GHz for J0741+312 
•^With a measured flux density of 2.67(±0. 10) Jy for J2052+365 
"■Wifli a measured flux density of 2.62(±0.07) Jy for J2052+365 

range, in our case at 6030 and 6035 MHz. With the con- 
figuration in A-array, three EVLA antennas were spread near- 
homogeneously over the North and West arm each. Though 
the East arm had a similar distribution of EVLA anteimas, the 
outer two were not operational at 6.0 GHz. 

For AU Gem, two hours were devoted to the 4750 MHz line 
followed by two hours on the 1665 and 1667 MHz lines of 
OH. Due to the lack of observations of an absolute standard 
calibrator (e.g., 3C286), the source 0741+312 was used for 
bandpass, phase and amplitude, and primary flux calibrations, 
using assumed flux densities of 2.0 Jy at 1.7 GHz and 1.6 Jy 
at 4.8 GHz as taken from the VLA calibrator list. System 
flux monitoring suggests this is in error with less than 20%. 
As a temporary inconvenience, onUne Doppler tracking was 
not used due to differences in the VLA and EVLA antenna 
control systems. Instead, the sky frequency was calculated 
from the LSR velocity using the NRAO online Dopset tool 
and the observations were taken in fixed-frequency mode. The 
LSR velocity scale is set to the mean sky frequency of the 
observations. Maximum deviations from this frequency (at 
the beginning and end of observing in each band) correspond 
to 0.50 channel widths at 4.8 GHz (< 0.1 km s"') and 0.36 
channel widths at 1.7 GHz (< 0.1 km s"'). 

Two 1.5 hour blocks were devoted to observations of the 
6030 and 6035 MHz lines in NML Cyg. By using only the 
EVLA antennas, online Doppler tracking was available. The 
data were calibrated in the standard way for VLA antennas us- 
ing 3C48 as standard flux calibrator and J2052+365 as phase 
cahbrator. The new EVLA system required more careful flag- 
ging of bad data than usual with the estabUshed VLA system. 
Data from some antennas were discarded entirely: one an- 
tenna did not have phase stability, while another antenna had 
no signal in the 6030 MHz LCP IF. A third antenna showed 
occasional phase jumps of about 60°. The absolute flux cal- 
ibration at 6.0 GHz has not yet been determined; flux den- 
sities and RMS noise levels given in Table 1 assume a scal- 
ing with frequency for 3C48 as calculated in the AIPS task 
SETJY (Greisen 2003). 

3. RESULTS 
3.1. AUGem 

Figure 1 shows the emission at 1665 and 1667 MHz in AU 
Gem. We detect ^ 250 mJy emission in each of the 1665 
and 1667 MHz transitions in the 12-13 kms"' range, corre- 
sponding to similarly bright 1667 MHz emission detected by 
Nguyen-Q-Rieu et al. (1979). Right circular polarized (RCP) 
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Fig. 2. — Manning smoothed spectra of tentative 6035 MHz emission 
vi'ithin 1" of the position of NML Cyg. Arbitrarily scaled 1612 MHz emis- 
sion from VLA archival data is shown with a dotted line for reference. 

4. DISCUSSION 

Theoretical work on OH pumping in model, well-behaved 
circumstellar shells predict that excited-state OH masers 
should not be seen, and observations in most circumstellar 
shells to date agree. Our non-detections of excited-state OH 
in AU Gem in the 4750 MHz line and in NML Cyg in the 
6030 and 6035 MHz lines support the prediction that excited- 
state OH should not be observable in circumstellar shells as 
part of their normal structure and evolution. Excepting tem- 
porary events that cannot easily be reobserved, the previous 
reports of detection of the excited-state OH lines in AU Gem 
and NML Cyg cannot be explained. 

4.1. AUGem 

Dismissing the not very sensitive observation by Jewell et 
al. (1985), our 10 mJy beam"' 3cr upper limit non-detection 
of the 4750 MHz excited-state OH emission contrasts the sole 
report of a ~100 mJy Su positive detection in this line by 
Claussen & Fix (1981). If not an unrelated source in the line- 
of-sight of the 0.'5 beam of the Arecibo telescope, this sug- 
gests two possibilities. 

The first is that the 4750 MHz maser emission in AU Gem 
has weakened substantially (by more than a factor of 10) be- 
tween 1981 and 2007. The lack of other known evolved stars 
with 4750 MHz OH masers precludes us from commenting 
on the phenomenology in this transition, but it is worth noting 
that in a comparable time period the 6035 MHz Une in NML 
Cyg weakened by more than a factor of 100 (Zuckerman et 
al. 1972; Desmurs et al. 2002, and discussed below). The sec- 
ond possibility is that the Claussen & Fix (1981) detection is 
spurious. It is based on a 5 ct spike in a single 0.6 km s"' spec- 
tral channel in frequency switched Arecibo data. It is possible 
that this emission arises from terrestrial interference or a bad 
correlator channel rather than a celestial source. The LSR 
velocity of this feature, 3.5 kms"', is outside the range of the 
dominant emission in either ground-state main line. However, 
it would still fall within the derived velocity range of the outer 
1667 MHz OH-peaks as argued in the first (Vexp = 14-15 km 
s"') case in Sect. 3.1. Even though the 1612 MHz OH, H2O 
and SiO transitions have not been detected (Fix & Weisberg 
1978; Nguyen-Q-Rieu et al. 1979; Olnon et al. 1980; Nyman 
et al. 1986), making AU Gem atypical for an (Type II) OH/IR 
star, we could not argue a special case scenario for AU Gem 
as for NML Cyg outlined below. AU Gem probably just does 
not have a circumstellar environment thick or dusty enough to 
sustain 1612 MHz maser emission and thus we think that it is 
unlikely to have ever had excited-state OH emission. 

4.2. NML Cyg 

The new non-detections of 6030 and 6035 MHz emission 
toward NML Cyg suggest two possibilities. Like AU Gem, 
this phenomenon could be a time-variable or single event 
in the history of NML Cyg as perhaps typical for any other 
OH/IR object during their evolution. However, as stated be- 



fore, because other searches for excited-state OH emission in 
circumstellar environments of a variety of evolved stars have 
yielded no detections, it appears that NML Cyg must some- 
how be special to have had 6.0 GHz excited-state OH maser 
emission. We reject the possibility of the original multi- 
channel ^ lOcr detection by Zuckerman et al. (1972) as being 
spurious. The second possible explanation is that NML Cyg 
indeed is special compared to other evolved OH/IR type stars 
due to the interaction with its environment. 

Since the discovery of emission in the ground-state OH 
lines in the optically obscured infrared star NML Cyg by Wil- 
son & Barrett (1968a,b), it has been extensively studied in 
the radio mainly in the very bright 1612 MHz satellite line 
of OH. The main lines at 1665 and 1667 MHz are present as 
well (Wilson & Barrett 1968a,b, and from VLA archive data), 
but with typically a factor of 100 less flux than the 1612 MHz 
emission, whereas the 1720 MHz satellite line was reported to 
be undetected by Wilson & Barrett (1968a,b). For compari- 
son in the following discussion we refer to the high resolution 
(MERLIN) ground-state blueshifted 1665 MHz mainline OH 
and fuU extent 1612 MHz OH data on NML Cyg described by 
Diamond et al. (1984) and Etoka & Diamond (2004). 

Though the asymmetric morphology in the 1612 MHz OH 
transition was modeled at first as a rotating disk at a position 
angle of about 150° toward the northwest (e.g., Masheder et 
al. 1974; Benson & Mutel 1979), subsequent interpretations 
favor a double 1612 MHz OH circumstellar shell (Diamond 
et al. 1984; Etoka & Diamond 2004) with the inner shell as 
a spherical expanding shell typical for an OH/IR star. The 
outer 1612 MHz OH shell only manifests itself as a curved 
arc-like shell-segment located to the northwest at about 2"3 
from the star at a velocity close to the stellar velocity, where 
the shell motion is predominantly tangential in the plane of the 
sky (Etoka & Diamond 2004, their Figures 8 and 9). The in- 
complete (blueshifted emission only) 1665 MHz mainline OH 
image shows a similar picture, with a less prominent arc lo- 
cated on the sky at about 1 ."3 to the northwest and between the 
two 1612 MHz shells (Etoka & Diamond 2004, their Figure 
7). It is noteworthy that the mass outflow is not constant; e.g., 
Danchi et al. (2001) deduce a 3.86 mas yr"' infrared proper 
motion of another double shell closer (< 0."3, versus 2."3 for 
the 1612 MHz OH arc) to the star, and an age difference be- 
tween the two inner expanding shells of 65 (±14) years. 

Toward the west-northwest, NML Cyg is surrounded by an 
H II region (Habing et al. 1982), originating from the ion- 
ization of outflowing material from NML Cyg by the intense 
photoionizing radiation field of the nearby Cyg OB2 associa- 
tion (Morris & Jura 1983). It is in particular interesting to note 
that Habing et al. (1982) overlay their H II observations on the 
red Palomar Sky Atlas E plate and identify a near-Unear fea- 
ture at about 30-35" west-northwest from NML Cyg. Habing 
et al. (1982) suggest that this depicts Ha emission, which we 
in turn recognize as a tracer for shocked material. Recent HST 
observations by Schuster et al. (2006) of the dust immediately 
surrounding NML Cyg outUnes the dissociation surface gen- 
erated by this radiation field, oriented in the same direction as 
the H n region and flie 1612 MHz OH arc (but with 0"25 ex- 
tent at a much smaller scale). We find it plausible that a bow 
shock-like front causes the OH molecules and dust toward the 
direction of the Cyg OB2 association to pile up, building up 
sufficient OH (colunm) density. The increased absorption of 
stellar radiation by the enhanced density of dust will radia- 
tively pump the 1612 MHz maser, causing a partial sheU or 
arc to appear at a line-of-sight velocity near the steUar veloc- 
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ity (Etoka & Diamond 2004). 

We note that 1720 MHz satellite OH masers generally are 
seen in high density shocks, as is Ha emission. However, re- 
cent calculations (Wardle 2007) show that excited-state OH 
might be observable at even higher densities. The variation 
in mass loss of NML Cyg may cause the northwest side pile 
up to occasionally be shocked and temporarily have increased 
density due to shells impacting on the near-stationary mate- 
rial. It could have been that Zuckerman et al. (1972) observed 
NML Cyg during such an impact event. Since then the shock 
likely has dissipated and the densities have fallen. If the 65 
years between the shells found by Danchi et al. (2001) is typi- 
cal, currently we are observing halfway between the impact of 
two shells, predicting another impact around the years 2030- 
2040. We speculate that while the excited-state emission is 
now gone, ground-state 1720 MHz OH emission may be de- 
tectable if the pumping is predominantly collisional. Obser- 
vations of the 1720 MHz transition at greater sensitivity (<C 
1 Jy) than that obtained by Wilson & Barrett (1968a,b) are 
required to test this hypothesis. 

Future discussions of maser pumping models should there- 
fore not be distracted by these two cases of reported excited- 
state emission in circumstellar environments of late-type stars 
anymore. Instead, current models seem sufficient to explain 
the pumping of masers in circumstellar shells of evolved stars 
as part of their normal structure and evolution. 

5. SUMMARY 

Although a temporary event in AU Gem causing 4750 MHz 
excited-state OH emission cannot be excluded, we conclude 



that it is probable that the original report by Claussen & Fix 
(1981) regards an unfortunate spurious detection and is not 
4750 MHz excited-state OH maser emission in the circum- 
stellar envirormient of AU Gem. 

The reported 6035 MHz excited-state OH emission by 
Zuckerman et al. (1972) in NML Cyg cannot be spurious and 
can be explained as originating from the special temporary 
conditions arising from a shock between the high mass out- 
flow of NML Cyg and the intense ionizing UV-radiation field 
of the nearby Cyg OB2 association. Although it is clear that 
the 6035 MHz emission in NML Cyg has weakened substan- 
tially, it is unclear whether the excited-state OH emission is 
still present at very low levels due to a lack of sensitivity in 
this early-stage of the VLA to EVLA transition. We suggest 
that the excited- state OH transitions in NML Cyg be reob- 
served with higher sensitivity when there are more EVLA an- 
tennas with this capability. We also recommend that the 1720 
MHz OH transition be observed as a tracer of a possible shock 
in the environment of NML Cyg. 

The transition from the VLA to the EVLA is well underway 
and already now offers great opportunities to observe at fre- 
quencies previously unavailable to the VLA user community. 

The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 
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